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Abstract 

In a previous paper (Gavignaud et al. 2006), we presented the type-1 Active Galactic Nuclei (AGN) sample obtained from the first epoch data 
of the VIMOS-VLT Deep Survey (VVDS). The sample consists of 130 faint, broad-line AGN with redshift up to z = 5 and 17.5 < I AB < 24.0, 
selected on the basis of their spectra. 

In this paper we present the measurement of the Optical Luminosity Function up to z = 3.6 derived from this sample, we compare our results with 
previous results from brighter samples both at low and at high redshift. 

Our data, more than one magnitude fainter than previous optical surveys, allow us to constrain the faint part of the luminosity function up to high 
redshift. 

By combining our faint VVDS sample with the large sample of bright AGN extracted from the SDSS DR3 (Richards et al., 2006b), we find that the 
model which better represents the combined luminosity functions, over a wide range of redshift and luminosity, is a luminosity dependent density 
evolution (LDDE) model, similar to those derived from the major X-surveys. Such a parameterization allows the redshift of the AGN space density 
peak to change as a function of luminosity and explains the excess of faint AGN that we find at 1.0 < z < 1.5. On the basis of this model we find, 
for the first time from the analysis of optically selected samples, that the peak of the AGN space density shifts significantly towards lower redshift 
going to lower luminosity objects. This result, already found in a number of X-ray selected samples of AGN, is consistent with a scenario of 'AGN 
cosmic downsizing", in which the density of more luminous AGN, possibly associated to more massive black holes, peaks earlier in the history of 
the Universe, than that of low luminosity ones. 

Key words, surveys-galaxies: high-redshift - AGN: luminosity function 



1. Introduction 

Active Galactic Nuclei (AGN) are relatively rare objects that ex- 
hibit some of the most extreme physical conditions and activity 
known in the universe. 

A useful way to statistically describe the AGN activity along 
the cosmic time is through the study of their luminosity func- 
tion, whose shape, normalization and evolution can be used to 
derive constraints on models of cosmological evolution of black 
holes (BH). At z<2.5, the luminosity function of optically se- 
lected type-1 AGN has been well studied since many years 
dBovle et all \l9M. iHewetFetail. [19911 IPel [1991 iBovle et al.L 
l2000t ICroom et all 120041) . It is usually described as a double 
power law, characterized by the evolutionary parameters L*(z) 
and 0*(z), which allow to distinguish between simple evolution- 
ary models such as Pure Luminosity Evolution (PLE) and Pure 
Density Evolution (PDE). Although the PLE and PDE mod- 
els should be mainly considered as mathematical descriptions 
of the evolution of the luminosity function, two different phys- 
ical interpretations can be associated to them: either a small 



fraction of bright galaxies harbor AGN, and the luminosities of 
these sources change systematically with time ('luminosity evo- 
lution'), or all bright galaxies harbor AGN, but at any given time 
most of them are in 'inactive' states. In the latter case, the frac- 
tion of galaxies with AGN in an 'active' state changes with time 
('density evolution'). Up to now, the PLE model is the preferred 
description for the evolution of optically selected QSOs, at least 

at low redshift (z < 2). 

W orks on high redshift type-1 AGN samples dWarren et al 
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IWolf et all I2003t iHunt et all |2004|) have shown that the 
number density of QSOs declines rapidly from z ~ 3 to z ~ 5. 
Since the size of complete and well studied samples of QSOs 
at high redshift is still relatively small, the rate of this decline 
and the shape of the high redshift luminosity function is not yet 
as we ll constrained as at low redshift. For example, iFan et al] 
d2001l) . studying a sample of 39 luminous high redshift QSOs at 
3.6 < z < 5.0, selected from the commissioning data of the Sloan 
Digital Sky Survey (SDSS), found that the slope of the bright end 
of the QSO luminosity function evolves with redshift, becoming 
flatter at high redshift, and that the QSO evolution from z = 2 
to z = 5 cannot be described as a pure luminosity evolution. A 
similar result on the flattening at high redshift of the slope of 
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the luminos ity function fo r lumin ous QSOs has been recently 
obtained bv lRichards et al.l (l2006bl) fr om the analysis of a m uch 
larger sample of SDSS QSOs (but see lFontanot et all d2007l) for 
different conclusions drawn on the basis of combined analysis of 
GOODS and SDSS QSOs). 

At the same time, a growing number of observations at differ- 
ent redshifts, in radio, optical and soft and hard X-ray bands, are 
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Hasinger et alll2005l). This evolution, now dubbed as "AGN cos 





mic downsizing" is described either as a direct evolution in the 
faint end slope or as "luminosity dependent density evolution" 
(LDDE), and it has been the subject of many speculations since 
it implies that the space density of low luminosity AGNs peaks 
at lower redshift than that of bright ones. 

It has been observed that, in addition to the well known 
local scale relations between the blac k hole (BH) masses and 
the p r operties of their host galaxies dKormendv & Richstonel 
119951; iMagorrian et all ll998llFerrarese & Merrittl |2000|) , also 
the galaxy spheroi d population fo l lows a similar pattern of "cos- 
mic downsizing" dCimatti et all 120061) . Various models have 
been proposed to explain this common evolutionary trend in 
AGN and spheroid galaxies. The majority of them propose that 
the feedback from the black hole growth plays a key r ole in 
determining the BH-ho st galaxy relations dSilk & ReesL Il998t 
iDi MatteoetaU 120051) and the co-evolution of black holes 
and their host galaxies. Indeed, AGN feedback can shut down 
the growt h of the most massive sys tems steepening the bright 
end slope CS cannapieco & OR|2004|) . while the feedback-driven 
QSO decay determine s the shape of the faint end of the QSO LF 
dHopkins et all 120061) . 

This evolutionary trend has not been clearly seen yet with 
optically selected type-1 AGN samples. By combining results 
from low and high redshifts, it is clear from the studies of op- 
tically selected samples that the cosmic QSO evolution shows a 
strong increase of the activity from z ~ out to z ~ 2, reaches a 
maximum around z - 2 - 3 and then declines, but the shape of 
the turnover and the redshift evolution of the peak in activity as 
a function of luminosity is still unclear. 

Most of the optically selected type-1 AGN samples stud- 
ied so far are obtained through various color selections of 
candidates, follow ed by spe ctroscopic confirmat ion (e.g. 2dF, 
ICroom etal.l2004l and SDSS. lRichards et al.l200"2l) , or grism and 
slitless spectroscopic surveys. These samples are expected to be 
highly complete, at least for luminous type-1 AGN, at either 
z < 2.2 or z > 3.6, where type-1 AGN show conspicuous colors 
in broad band color searches, but less complete in the redshift 
range 2.2 < z < 3.6 (Richards et al. 2002). 

An improvement in the multi-color selection in optical bands 
is through the simultaneous use of m any broad and me dium band 
filters as in the COMBO-17 survey dWolf et alll2003l This sur- 
vey is the only optical survey so far which, in addition to cov- 
ering a redshift range large enough to see the peak of AGN ac- 
tivity, is also deep enough to sample up to high redshift type-1 
AGN with luminosity below the break in the luminosity func- 
tion. However, only photometric redshifts are available for this 
sample and, because of their selection criteria, it is incomplete 
for objects with a small ratio between the nuclear flux and the 
total host galaxy flux and for AGN with anomalous colors, such 
as, for example, the broad absorption line (BAL) QSOs , which 
have on average redder c olors and account for ~ 10 - 15 % of the 
overall AGN population dHewett & Foltzl 120031) . 



The VIMOS-VLT Deep Survey dLe Fevre et all 120051) is a 
spectroscopic survey in which the target selection is purely flux 
limited (in the I-band), with no additional selection criterion. 
This allows the selection of a spectroscopic type-1 AGN sample 
free of color and/or morphological biases in the redshift range 
z > 1. An obvious advantage of such a selection is the possi- 
bility to test the compl eteness of the most current surveys (see 
iGavignaud et allr2006l Paper I), based on morphological and/or 
color pre- selection, and to study the evolution of type-1 AGN 
activity in a large redshift range. 

In this paper we use the type-1 AGN sample selected from 
the VVDS to derive the luminosity function in the redshift range 
1 < z < 3.6. The VVDS type-1 AGN sample is more than one 
magnitude deeper than any previous optically selected sample 
and allow thus to explore the faint part of the luminosity func- 
tion. Moreover, by combining this LF with measurement of the 
LF in much larger, but very shallow, surveys, we find an analyt- 
ical form to dercribe, in a large luminosity range, the evolution 
of type-1 AGN in the redshift range 0< z <4. The paper is or- 
ganized as follows: in Section 2 and 3 we describe the sample 
and its color properties. In Section 4 we present the method used 
to derive the luminosity function, while in Section 5 we com- 
pare it with previous works both at low and high redshifts. The 
bolometric LF and the comparison with the results derived from 
samples selected in different bands (from X-ray to IR) is then 
presented in Section 6. The derived LF fitting models are pre- 
sented in Section 7 while the AGN activity as a function of red- 
shift is shown in Section 8. Finally in section 9 we summarize 
our results. Throughout this paper, unless stated otherwise, we 
assume a cosmology with Q. m = 0.3, Qa = 0.7 and Ho = 70 km 
s _1 Mpc~'. 



2. The sample 

Our AGN sample is extracted from the first ep och data of the 
VIMO S-VLT Deep Survey, performed in 2002 dLe Fevre et all 
120051) . 

The VVDS is a spectroscopic survey designed to measure 
about 150,000 redshifts of galaxies, in the redshift range < z < 
5, selected, nearly randomly, from an imaging survey (which 
consists of observations in U, B, V, R and I bands and, in a 
small area, also K-band) designed for this pu rpose. Full de- 
tails ab o ut VIMOS photometry ca n be found in iLe Fevre et al.l 
(l2004al) . iMcCracken et alJ d2003l). iRadovich et alj d2004l) for 
the U-band and llovino etaU (120051) for the K -band. In this 
work we will as well use the Gal ex UV-catalog dArnouts et all 
l2005t ISchiminovich et all 120051) . the u*,g',r',i',z' photometry 
obtained in the frame of the Canada-France-Hawai i Legacy 
Survey (CFHTLSfl UKIDSS dLawrence et all [20061) . and the 
Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) 
dLonsdale et all 120031 12004|) . The spectroscopic VVDS survey 
consists of a deep and a wide survey and it is based on a sim- 
ple selection function. The sample is selected only on the basis 
of the I band magnitude: 17.5 < Iab < 22.5 for the wide and 
17.5 < Iab < 24.0 for the deep sample. For a detailed descrip- 
tio n of the spectroscopic survey strategy and the first epoch data 
see lLe Fevre et all d2005l) . 

Our sample consists of 130 AGN with < z < 5, selected 
in 3 VVDS fields (0226-04, 1003+ 01 and 2217-00) and i n the 
Chandra Deep Field South (CDFS. EeFevre et all l2004bl) . All 
of them are selected as AGN only on the basis of their spectra, 
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Figure 1. Distribution of absolute magnitudes and redshifts of 
the total AGN sample. Open circles are the objects with am- 
biguous redshift, shown at all their possible z values. The dotted 
and dashed lines represent the magnitude limits of the samples: 
Iab < 22.5 for the wide sample and Iab < 24.0 for the deep 
sample. 
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Figure 2. Composite spectra derived for our AGN with se- 
cure redshift in the 02h field, divided in a "bright" (19 objects 
at M1450 <-22.15, dotted curve) and a "faint" (31 objects at 
M1450 >-22.15, dashed curve) sample. We consider here only 
AGN with z > 1 (i.e. the AGN used in to compute the lumi- 
nosity function). The SDSS composite spectrum is shown with 
a solid line for comparison. 



irrespective of their morphological or color properties. In partic- 
ular, we selected them on the basis of the presence of at least one 
broad emission line. We discovered 74 of them in the deep fields 
(62 in the 02h field and 12 in the CDFS) and 56 in the wide fields 
(18 in the lOh field and 38 in the 22h field). This represents an 
unprecedented complete sample of faint AGN, free of morpho- 
logical or color selection bias. The spectroscopic area covered 
by the First Epoch Data is 0.62 deg 2 in the deep fields (02h field 
and CDFS) and 1.1 deg 2 in the wide fields (lOh and 22h fields). 

To each object we have assigned a value for the spectro- 
scopic redshift and a spectroscopic quality flag which quantifies 
our confidence level in that given redshift. As of today, we have 
115 AGN with secure redshift, and 15 AGN with two or more 
possible values for the redshift. For these objects, we have two 
or more possible redshifts because only one broad emission line, 
with no other narrow lines and/or additional features, is detected 
in the spectral wavelength range adopted in the VVDS (5500 - 
9500 A) (see Figure 1 in Paper I). For all of them, however, a 
best solution is proposed. In the original VVDS AGN sample, 
the number of AGN with this redshift degeneracy was 42. To 
solve this problem, we have first looked for the objects already 
observed in other spectroscopic surveys in the same areas, solv- 
ing the redshift for 3 of them. For the remaining objetcs, we 
performed a spectroscopic follow-up with FORS1 on the VLT 
Unit Telescope 2 (UT2). With these additional observations we 
found a secure redshift for 24 of our AGN with ambiguous red- 
shift determination and, moreover, we found that our proposed 
best solution was the correct one in ~ 80% of the cases. On the 
basis of this result, we decided to use, in the following analysis, 
our best estimate of the redshift for the small remaining fraction 
of AGN with ambiguous redshift determination (15 AGN). 

In Figure Q] we show the absolute B-magnitude and the 
redshift distributions of the sample. As shown in this Figure, 
our sample spans a large range of luminosities and consists of 
both Seyfert galaxies (M B >-23; -59%) and QSOs (M B <-23; 
~41%). A more detailed and exhaustive description of the prop- 



erties of the AGN sample is given in Paper I dGavignaud et all 
120061) and the complete list of BLAGN in our wide and deep 
samples is available a s an electronic Table in Appendix of 
iGavignaud et alj d2006l) . 

3. Colors of BLAGNs 

As already discussed in Paper I, the VVDS AGN sample shows, 
on average, redder colors than those expected by comparing 
them, for example, w ith the color track derived from the SDSS 
composite spectrum dVanden Berk et all 1200 lb . In Paper I we 
proposed three possible explanations: (a) the contamination of 
the host galaxy is reddening the observed colors of faint AGN; 
(b) BLAGN are intrinsically redder when they are faint; (c) the 
reddest colors are due to dust extinction. On the basis of the sta- 
tistical properties of the sample, we concluded that hypothesis 

(a) was likely to be the more correct, as expected from the faint 
absolute magnitudes sampled by our survey, even if hypotheses 

(b) and (c) could not be ruled out. 

In Figure|2]we show the composite spectra derived from the 
sample of AGN with secure redshift in the 02h field, divided 
in a "bright" and a "faint" sample at the absolute magnitude 
A^i450 = -22.15. We consider here only AGN with z > 1, which 
correspond to the AGN used in Section [4] to compute the lumi- 
nosity function. The choice of the reference wavelength for the 
absolute magnitude, A = 1450 A, is motivated by our photo- 
metric coverage. In fact, for most of the objects it is possible 
to interpolate M1450 directly from the observed magnitudes. In 
the same plot we show also the SDSS composite spectrum (solid 
curve) for comparison. Even if also the "bright" VVDS compos- 
ite (dotted curve) is somewhat redder than the SDSS one, it is 
clear from this plot that the main differences occur for faintest 
objects (dashed curve). 

A similar result is shown for the same sample in the upper 
panel of Figure [3] where we plot the spectral index a as a func- 
tion of the AGN luminosity. The spectral index is derived here 
by fitting a simple power law /(v) = v~ a to our photometric data 
points. This analysis has been performed only on the 02h deep 
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Figure 3. Upper Panel: Distribution of the spectral index a as 
a function of M1450 for the same sample of AGN as in Figure 
[2] The spectral index is derived here by fitting a simple power 
law /(v) = v~ a to our photometric data points. Asterisks are 
AGN morphologically classified as extended and the grey point 
is a BAL AGN. Bottom Panels: Distribution of the spectral in- 
dex a for the same sample of AGN. All the AGN in this sample 
are shown in the first of the three panels, while the AGN in the 
"bright" and "faint" sub-samples are shown in the second and 
third panel, respectively. The dotted curve in the second panel 
corresponds to the gaussian fit of the bright sub-sample and it is 
reported also in the third panel to highlight the differences in the 
a distributions of the two sub-samples. 



sample, since for the wide sample we do not have enough photo- 
metric coverage to reliably derive the spectral index. Most of the 
AGN with a > 1 are fainter than M1450 = -22.15, showing that, 
indeed, the faintest objects have on average redder colors than 
the brightest ones. The outlier (the brightest object with large a, 
i.e. very red colors, in the upper right corner of the plot) is a BAL 
AGN. 

The three bottom panels of Figure [3] show the histograms of 
the resulting power law slopes for the same AGN sample. The 
total sample is plotted in the first panel, while the bright and the 
faint sub-samples are plotted in the second and third panels, re- 
spectively. A Gaussian curve with < a >= 0.94 and dispersion 
cr = 0.38 is a good representation for the distribution of about 
80% (40/50) of the objects in the first panel. In addition, there 
is a significant tail (~ 20%) of redder AGN with slopes in the 
range from 1.8 up to ~ 3.0. The average slope of the total sample 
(~ 0.94) is redder than the fit to the SDSS composite (~ 0.44). 
Moreover, the distribution of a is shifted toward much larger val- 



ues (red der continua) than th e similar distribution in the SDSS 
sample dRichards et all [2003). For example, only 6% of the ob- 
jects in the SDSS sample have a > 1.0, while this percentage is 
57% in our sample. 

The differences with respect to the SDSS sample can be 
partly due to the differences in absolute magn itude of the two 
samp les (M, <-22.0 for the SDSS sample dSchneider et all 
2003) and M B <-20.0 for the VVDS sample). In fact, if we con- 
sider the VVDS "bright" sub-sample, the average spectral index 
< a > becomes ~ 0.71, which is closer to the SDSS value (even 
if it is still somewhat redder), and only two objects (~8% of the 
sample) show values not consistent with a gaussian distribution 
with cr ~0.32. Moreover, only 30% of this sample have a > 1.0. 

Most of the bright SDSS AGNs with a > 1 are interpreted by 
iRichards et alj d2003l) to be dust-reddened, althoug h a fraction 
of the m is likely to be due to intrinsically red AGN dHall et all 
120061) . At fainter magnitude one would expect both a larger frac- 
tion of dust-re ddened objects (in analogy with indications from 
the X-ray data dBrandt et all l2000t iMushotzkv et all l2000h and 
a more significant contamination from the host galaxy. 

We have tested these possibilities by examining the global 
Spectral Energy Distribution (SED) of each object and fitting 
the observed fluxes f t, s with a combination of AGN and galaxy 
emission, allowing also for the possibility of extinction of the 
AGN flux. Thanks to the multi-wavelength coverage in the deep 
field in which we have, in addition to VVDS bands, also data 
from GALEX, CFHTLS, UKIDSS and SWIRE, we can study 
the spectral energy distribution of the single objects. In particu- 
lar, we assume that: 



fobs = C.Ucn ■ W-°+ A * + C 2 f CA L (1) 



and, using a library of galaxy and AGN templates, we find the 
best parameters c\, ci and Eg_y for e a ch obje ct. We used the 
AGN SED derived by IRichards et al.l d2006al) with an SMC- 
like dust-reddening law dPrevot et all 1 19841) with the form 
A^/Eb-v = 1.39/ t ^j; 2 , a nd a library of galaxies template by 
iBruzual & Charlotl d2003l) . 

We found that for ~37% of the objects, the observed flux 
is fitted by a typical AGN power law (pure AGN), while 44% 
of the sources require the presence of a contribution from the 
host galaxy to reproduce the observed flux. Only 4% of the ob- 
jects are fitted by pure AGN + dust, while the remaining 15% 
of objects require instead both contributions (host galaxy con- 
tamination and presence of dust). As expected, if we restrict the 
analysis to the bright sample, the percentage of pure AGN in- 
creases to 68%, with the rest of the objects requiring either some 
contribution from the host galaxy (~2 1 %) or the presence of dust 
oscuration (~ 1 1 %). 

In Figure [4] we show 4 examples of the resulting fits: (i) pure 
AGN; (ii) dust-extincted AGN; (iii) AGN contaminated by the 
host galaxy; (iv) dust-extincted AGN and contaminated by the 
host galaxy. The dotted line corresponds to the AGN template 
before applying the extinction law, while the solid blue line cor- 
responds to the same template, but extincted for the given Eb-v', 
the red line corresponds to the galaxy template and, finally, the 
black line is the resulting best fit to the SED. The host galaxy 
contaminations will be taken into account in the computation of 
the AGN absolute magnitude for the luminosity function. 
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Figure 4. Four examples of different decompositions of the ob- 
servedSEDs of our objects. Since for A < 1216A, corresponding 
to the Lya line, the observed flux is expected to decrease because 
of intervening absorption, all the photometric data at A < 12 16 A 
are not considered in the fitting. The only requested constraint is 
that they lie below the fit. The four fits shown in this Figure cor- 
respond, from top to bottom, to pure- AGN, dust-extincted AGN, 
AGN and host galaxy, dust-extincted AGN and host galaxy. The 
dotted line corresponds to the AGN template before applying 
the extinction law, while the solid blue line corresponds to the 
same template, but extincted for the given Eg-v- The red line 
(third and fourth panel) corresponds to the galaxy template and, 
finally, the black line is the resulting best fit to the SED. Arrows 
correspond to 5<x upper limits in case of non detection in the IR. 



4. Luminosity function 

4. ) . Definition of the redshift range 

For the study of the LF we decided to exclude AGN with z < 
1.0. This choice is due to the fact that for 0.5 < z < 1.0 the 
only visible broad line in the VVDS spectra is H/3 (see Figure 
1 of Paper I). This means that all objects with narro w or almost 
narro w H/3 and broad Ha (type 1.8, 1.9 AGN; see lOsterbrockl 
1 198 lb would not be included in our sample, because we include 
in the AGN sample all the objects with at least one visible broad 
line. Since at low luminosities the number of intermediate type 
AGN is not negligible, this redshift bin is likely to be under- 
populated and the results would not be meaningful. 

At z < 0.5, in principle we have less problems, because also 
Ha is within the wavelength range of the VVDS spectra, but, 
since at this low redshift, our sampled volume is relatively small 
and QSOs rare, only 3 objects have secure redshifts in this red- 
shift bin in the current sample. For these reasons, our luminosity 
function has been computed only for z > 1.0 AGN. As already 
mentioned in Section 2, the small fraction of objects with an am- 



biguous redshift determination have been included in the compu- 
tation of the luminosity function assuming that our best estimate 
of their redshift is correct. 

The resulting sample used in the computation of the LF consists 
thus of 121 objects at 1< z <4. 

4.2. Incompleteness function 

Our incompleteness function is made up of two terms linked, re- 
spectively, to the selection algorithm and to the spectral analysis: 
the Target Sampling Rate (T SR) and the Spectr oscopic Success 
Rate (SSR) defined following llibert et alj (120051) . 

The Target Sampling Rate, namely the ratio between the ob- 
served sources and the total number of objects in the photometric 
catalog, quantifies the incompleteness due to the adopted spec- 
troscopic selection criterion. The TSR is similar in the wide and 
deep sample and runs from 20% to 30%. 

The Spectroscopic Success Rate is the probability of a spec- 
troscopically targeted object to be securely identified. It is a com- 
plex function of the BLAGN redshift, apparent magnitude and 
intrinsic spectral energy distribution and it has been estimated 
by simulating 20 Vimos pointings, for a total of 2745 spectra. 

Full details on TS R and SSR can be found in Paper I 
dGavignaud et aUl2006l) . We account for them by computing for 
each object the associated weights w tsr = 1 /TSR and w ssr = 
1/SSR; the total weighted contribution of each object to the 
luminosity function is then the product of the derived weights 
(w' sr x w ssr ). 

4.3. Estimate of the absolute magnitude 

We derived the absolute magnitude in the reference band from 
the apparent magnitude in the observed band as: 

M = m obs - 5log w (dl(z)) -25 -k (2) 

where M is computed in the band in which we want to compute 
the luminosity function, m \, s is the observed band from which 
we want to calculate it, dl(z) is the luminosity distance expressed 
in Mpc and k is the k-correction in the reference band. To make 
easier the comparison with previous results in the literature, we 
computed the luminosity function in the B-band. 

To minimize the uncertainties in the adopted k-correction, 
m t, s for each object should be chosen in the observed band 
which is sampling the rest-wavelength closer to the band in 
which the luminosity function is computed. For our sample, 
which consists only of z > 1 objects, the best bands to use to 
compute the B-band absolute magnitudes should be respectively 
the I-, J- and K-bands going to higher redshift. Since however, 
the only observed band available for the entire sample (deep and 
wide), is the I-band, we decided to use it for all objects to com- 
pute the B-band magnitudes. This means that for z > 2, we 
introduce an uncertainty in the absolute magnitudes due to the 
k-correction. We computed the absolute ma gnitude considering 
the te mplate derived from the SDSS sample (IVanden Berk et all 
|200T1) . 

As discussed in Section [3] the VVDS AGN sample shows 
redder colors than those typical of normal, more luminous AGN 
and this can be due to the combination of the host galaxy contri- 
bution and the presence of dust. Since, in this redshift range, the 
fractional contribution from the host galaxies is expected to be 
more significant in the I-band than in bluer bands, the luminos- 
ity derived using the I-band observed magnitude could, in some 
cases, be somewhat overestimated due to the contribution of the 
host galaxy component. 
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Figure 5. Real (full circles; AGN in the deep sample) and simu- 
lated (open triangles; AGN in the wide sample) B-band absolute 
magnitude differences as a function of M B (TOT) (upper panel) 
and redshift (bottom panel). Mb(TOT) is the absolute magnitude 
computed considering the total observed flux, while M#(AGN) 
is the absolute magnitude computed after subtracting the host- 
galaxy contribution. 



We estimated the possible impact of this effect on our re- 
sults in the following way. From the results of the analysis of the 
SED of the single objects in the deep sample (see Sectional we 
computed for each object the difference mi(TOT) - mj(AGN) 
and, consequently, M B (TOT) - Mb(AGN). This could allow us 
to derive the LF using directly the derived Mb(AGN), resolv- 
ing the possible bias introduced by the host galaxy contami- 
nation. These differences are shown as full circles in Figure 
[5] as a function of absolute magnitude (upper panel) and red- 
shift (lower panel). For most of the objects the resulting dif- 
ferences between the total and the AGN magnitudes are small 
(AM<0.2). However, for a not negligible fraction of the faintest 
objects (Mg >-22.5, z <2.0) these differences can be signifi- 
cant (up to ~1 mag). For the wide sample, for which the more 
restricted photometric coverage does not allow a detailed SED 
analysis and decomposition, we used simulated differences to 
derive the Mb(AGN). These simulated differences have been de- 
rived through a Monte Carlo simulation on the basis of the bi- 
variate distribution AM(M,z) estimated from the objects in the 
deep sample. AM(M,z) takes into account the probability distri- 
bution of AM as a function of Mb and z, between and the solid 
line in Figure [5] derived as the envelope suggested by the black 
dots. The resulting simulated differences for the objects in the 
wide sample are shown as open triangles in the two panels of 
Figure [5] 

The AGN magnitudes and the limiting magnitudes of the 
samples have been corrected also for galactic extinction on the 
basis of the mean extinctio n values E(B-V) in each field derived 
from lSchlegel et al.l d!998l) . Only for the 22h field, where the ex- 
tinction is highly variable across the field, we used the extinction 
on the basis of the position of individual objects. The resulting 
corrections in the I-band magnitude are A/ a* 0.027 in the 2h and 
lOh fields and A/ = 0.0089 in the CDFS field, while the average 
value in the 22h field is A/ = 0.065. These corrections have been 
applied also to the limiting magnitude of each field. 



4.4. The 1/V max estimator 

We derived the binned representat ion of the lumin osity function 
using the usual 1 / V maA estimator (Sc hmidtl 1 19681) . which gives 
the space density contribution of individual objects. The lumi- 
nosity function, for each redshift bin (z - Az/2 ; z + Az/2), is 
then computed as: 



<d(m) = — y 

AM 



I M+AM/2 w tsr w ssr 



M-AM/2 



V 



(3) 



where Vmaxj is the comoving volume within which the ; ob- 
ject would still be included in the sample. w' sr and wj sr are re- 
spectively the inverse of the TSR and of the SSR, associated to 
the i' h object. The st atistical uncertainty on <t>(M) is given by 
iMarshall et ail (fl983l) : 



AM \ 



(w ts 



M+AM/2 
M-AM/2 max,i 



') 2 



(4) 



We combined o ur samples at d i fferent depths using the 
method proposed by lAvni & Bahcalll dl980b . In this method it 
is assumed that each object, characterized by an observed red- 
shift z, and intrinsic luminosity L,, could have been found in any 
of the survey areas for which its observed magnitude is brighter 
than the corresponding flux limit. This means that, for our total 
sample, we consider an area of: 

n tot (m) = a deep +Q wide = 1.72 deg 2 for 1 7 .5 < I AB < 22 .5 
and 

Q tof (m) = Q deep = 0.62 deg 2 for 22.5 < I AB < 24.0 

The resulting luminosity functions in different redshift 
ranges are plotted in Figure|6]and|2l where all bins which contain 
at least one object are plotted. The LF values, together with their 
lcr errors and the numbers of objects in each absolute magnitude 
bin are presented in Table Q] The values reported in Table Q] and 
plotted in Figures [6] and [7] are not corrected for the host galaxy 
contribution. We have in fact a posteriori verified that, even if the 
differences between the total absolute magnitudes and the mag- 
nitudes corrected for the host galaxy contribution (see Section 
14.3b can be significant for a fraction of the faintest objects, the 
resulting luminosity functions computed by using these two sets 
of absolute magnitudes are not significantly different. For this 
reason and for a more direct comparison with previous works, 
the results on the luminosity function presented in the next sec- 
tion are those obtained using the total magnitudes. 



5. Comparison with the results from other optical 
surveys 

We derived the luminosity function in the redshift range 1.0< 
z <3.6 and we compared it with the results from other surveys at 
both low and high redshift. 

5. 1 . The low redshift luminosity function 

In Figure [6] we present our luminosity function up to z = 2.1. 
The Figure show our LF data points (full circles) derived in two 
redshift bins: 1.0 < z < 1.55 and 1.55 < z < 2.1 c ompared with 
the LF fits de rived from the 2dF QSO sample by ICroom et al.l 
(12004 and bv lBovle et all d2000l) . with the COMBO- 17 sample 
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Figure 6. Ou r rest-frame B-band l uminosity function , derived in the reds hift bins 1 .0 < z < 1 .55 and 1 .55 < z < 2.1, compared with 
the 2dFQRS dCroom et all l2004t iBovle et all l2000l). COMBO-17 data dWolfet all [2003h and with the 2dF-SDSS (2SLAQ) data 
jRichards et all 120051) . The curves in the Figure show the PLE fit models derived by these authors. The thick parts of the curves 
correspond to the luminosity range covered by the data in each sample, while the thin parts are extrapolations based on the best fit 
parameters of the models. 



Wolf et al.1 d2003l). an d with the 2dF-SDSS (2SLAQ) LF fit 
Richards et al. In each panel the curves, computed 



by 
by 

for the average z of the redshift range, correspond to a double 
power law luminosity function in which the evolution with red- 
shift is characterized by a pure luminosity evolution modeled as 
M* b (z) = M*(0) - 2.5(/tiz + k 2 z 2 ). Moreover, the thick parts of the 
curves show the luminosity range covered by the data in each of 
the comparison samples, while the thin parts are extrapolation 
based on the the best fit parameters of the models. 

We start considering the comparison with the 2dF and the 
COMBO-17 LF fits. As shown in Figure[6] our bright LF data 
points connect rather smoothly to the faint part of the 2dF data. 
However, our sample is more than two magnitudes deeper than 
the 2dF sample. For this reason, a comparison at low luminosity 
is possible only with the extrapolations of the LF fit. At z > 1 -55, 
while the Boyle's model fits well our faint LF data points, the 
Croom's extrapolation, being very flat, tends to underestimate 
our low luminosity data points. At z < 1.55 the comparison is 
worse: as in the higher redshift bin, the Boyle's model fits our 
data better than the Croom's one but, in this redshift bin, our 
data points show an excess at low luminosity also with respect to 
Boyle's fit. This trend is similar to what shown also by the com- 
parison with the fit of the COMBO-17 data which, differently 
from the 2dF data, have a low luminosity limit closer to ours: at 
z > 1.55 the agreement is very good, but in the first redshift bin 
our data show again an excess at low luminosity. This excess is 
likely due to the fact that, because of its selection criteria, the 
COMBO-17 sample is expected to be significantly incomplete 
for objects in which the ratio between the nuclear flux and the 
total host galaxy flux is sma ll. Finally, we compare our data with 
the 2SLAQ fits derived by iRichards et all (120051) . The 2SLAQ 
data are derived from a sample of AGN selected from the SDSS, 



at 18.0 < g < 21.85 and z < 3, and observed with the 2-degree 
field instrument. Similarly to the 2dF sample, also for this sam- 
ple the LF is derived only for z < 2.1 and Mb < -22.5. The plot- 
ted dot-dashed curve corresponds to a PLE model in which they 
fix ed most of the para meters of the model at the values found 
bv lCroom et al.l d2004l) . leaving to vary only the faint end slope 
and the normalization constant cD*. In this case, the agreement 
with our data points at z < 1 .55 is very good also at low lu- 
minosity. The faint end slope f ound in this case is /3 = -1.45, 
which is similar to that found bv lBovle et ai1(l2000h (B = -1.58) 
and significantly steeper than that found bylCroom et al.l (120041) 
(J3 = -1.09). At z > 1.55, the IRichards et al.l(l2005l) LF fit tends 
to overestimate our data points at the faint end of the LF, which 
suggest a flatter slope in this redshift bin. 



The first conclusion from this comparison is that, at low red- 
shift (i.e. z < 2.1), the data from our sample, which is ~2 mag 
fainter than the previous spectroscopically confirmed samples, 
are not well fitted simultaneously in the two analyzed redshift 
bins by the PLE models derived from the previous samples. 
Qualitatively, the main reason for this appears to be the fact that 
our data suggest a change in the faint end slope of the LF, which 
appears to flatten with increasing redsh ift. This trend, already 
highlighted by previous X-ray surv eys (lLa Franca et all 120021: 
lUeda et all 120031 iFiore et all 120031) suggests that a simple PLE 
parameterization may not be a good representation of the evolu- 
tion of the AGN luminosity function over a wide range of red- 
shift and luminosity. Different model fits will be discussed in 
Section 7. 
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Figure 7. Our luminosity function, at 1450 A rest-frame, in the 
redshift ran ge 2.1<z<3 . 6, com pared with data from other high- 
z samples ( Hunt et al.l d2004l) at z — 3; Combo- 17 data from 
IWolf et alJ d2003l) at 2.4 < z < 3.6; data from IWarren et al 
d 19941) at 2.2 < z < 3.5 and the SDSS data from iFan et al 
(boOll) ). The SDSS data points at 3.6< z <3.9 have been 



e volved to z-3 usin g the luminosity evolution of iPeil d 19951) as 
in lHunt et all (|2004). The curves show some model fits in which 
the thick parts of the curves correspond to the luminosity range 
covered by the data samples, while the thin parts are model ex- 
trapolation. For this plot, an Q,„ = 1, Qa = 0, h = 0.5 cosmology 
has been assumed for comparison with the previous works. 



5.2. The high redshift luminosity function 

The comparison of our LF data points for 2.1< z <3.6 (full 
circles) with the results from other samples in similar redshift 
ranges is shown in Figure [7j In this Figure an Q.,,, = 1, Oa = 0, 
h = 0.5 cosmology has been assumed for comparison with pre- 
vious works, and the absolute magnitude has been computed at 
1450 A. As before, the thick parts of the curves show the lu- 
minosity ranges covered by the various data samples, while the 
thin parts are model extrapolations. In terms of number of ob- 
jects, depth and covered area, th e only sample com parable to 
ours is the COMBO- 17 sample dWolf etall 120031) . which, in 
this redshift range, consists of 60 AGN candidates over 0.78 
square degree. At a similar depth, in terms of abs olute mag- 
nitude , we show also the data from the sample of iHunt et al.l 
d2004l) . which however co nsists of 1 1 AGN in the redshift range 
< z > ±cr- =3.03+0.35 dSteidel et all 120021) . Given the small 
number of objects, the corresponding Hu nt model fit was de - 
rived including also the Warre n data point s dWarren et al.Lll994b . 
Moreover, they assumed the IPeil d 19951) luminosity evolution 
model, adopting the same values for L* and O*, leaving free to 
vary the two slopes, both at the faint and at the bright end of 
the LF. For comparison we show also the original Pei model 
fit derived from the empiric al lu minosity func t ion es timated 
by lHartwick & Schadd d!990l) and IWarren et all d!994l) . In the 



same plot we show also the model fit derived from a sample of 
~100 z ~ 3 (U-dropout) Q SO candidates by Siana et al. (pri- 
vate comunication; see also lSiana et al.ll2006l) . This sample has 
been selected by using a simple optical/lR photometric selec- 
tion at 19< r 1 <22 and the model fit has been derived by fix- 
ing the bright end slope at z=-2.85 as determined by SDSS data 
dRichardsetaUl2006bl) . 

In general, the comparison of the VVDS data points with 
those from the other surveys shown in Figure [7j shows a satis- 
factory agreement in the region of overlapping magnitudes. The 
best model fit which reproduce our LF data points at z ~ 3 is the 
Siana model with a faint end slope /? = -1.45. It is interesting to 
note that, in the faint part of the LF, our data points appear to be 
higher with respect to the IHunt et all d2004l) fit and are instead 
closer to the ext rapolation of the o riginal Pei model fit. This dif- 
ference with the lHunt et al.l (12004) fit is probably due to the fact 
that, having only 1 1 AGN in their faint sample, their best fit to 
the faint-end slope was poorly constrained. 



6. The bolometric luminosity function 

The comparison between the AGN LFs derived from samples 
selected in different bands has been for a long time a critical 
point in the studies o f the AGN luminosity function. Recently, 
iHopkins et al.l d2007l) . combining a large number of LF measure- 
ments obtained in different redshift ranges, observed wavelength 
bands and luminosity intervals, derived the Bolometric QSO 
Luminosity Function in the redshift range z = - 6. For each 
observational band, they derived appropriate bolometric correc- 
tions, taking into account the variation with luminosity of both 
the average absorption properties (e.g. the QSO column density 
Nh from X-ray data) and the average global spectral energy dis- 
tributions. They show that, with these bolometric corrections, it 
is possible to find a good agreement between results from all 
different sets of data. 

We applied to our LF d ata points the bolome tric corrections 
given by Eqs. (2) and (4) of lHopkins et all d2007l) for the B-band 
and we derived the bolometric LF shown as black dots in Figure 
[8] The s olid line represents the bolometric LF best fit model de- 
rived by IHopkins et al.l d2007l) and the colored data points cor- 
respond to different samples: green points are from optical LFs, 
blue and red points are from soft-X and hard-X LFs, respec- 
tively, and finally the cyan points are from the mid-IR LFs. All 
these bolometric LFs data poin ts have been derived following 
the same procedure described in IHopkins et all d2007l) . 

Our data, which sample the faint part of the bolometric lu- 
minosity function better than all previous optically selected sam- 
ples, are in good agreement with all the other samples, selected 
in different bands. Only in the last redshift bin, our data are quite 
higher with respect to the samples selected in other wavelength 
bands. The agreement remains however good with the COMBO- 
17 sample which is the only optically selected sample plotted 
here. This effect can be attributed to the fact that the conversions 
used to compute the Bolometric LF, being derived expecially for 
AGN at low redshifts, become less accurate at high redshift. 

Our data show moreover good agr eement also with the 
model fit derived bv lHopkins et al.l d2007l) . By try ing various an- 
alytic fits to the bolometric luminosity function IHopkins et al.l 
d2007l) concluded that neither pure luminosity nor pure density 
evolution represent well all the data. An improved fit can in- 
stead be obtained with a luminosity dependent density evolution 
model (LDDE) or, even better, with a PLE model in which both 
the bright- and the faint-end slopes evolve with redshift. Both 
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Figure 8. Bolometric luminosity function derived in three 
redshift bins from our data (black dots), compared with 
iHopkins et al ] (I2007h best-fit model and the data-sets used in 
their work. The central redshift of each bin is in dicated in each 
panel. Here, we adopted the same color-code as in lHopkins et al.l 
(120071) . but for more clarity we limited the number of samples 
presented in the Figure. Red symbol s corres pond to hard X-ra y 
surveys (squares: iBarger et al.ll2005t circles: lUeda et al.ll2.QQ3b . 



Blue to soft X-ray surveys (squares: ISilverman et al. 20051: cir- 
cles: Hasinger et all 120051) . Cyan to infra-red surveys (circles: 
iBrown et al.ll2006t squares: iMatute et al .1120061) . For the optical 
surveys we are showing here, with green circles, the data from 
the COMBO- 17 survey dWolf et all 120031) . which is comparable 
in depth to our sample. 



these models can reproduce the observed flattening with redshift 
of the faint end of the luminosity function. 



7. Model fitting 

In this Section we discuss the results of a number of different fits 
to our data as a function of luminosity and redshift. For this pur- 
pose, we computed the luminosity function in 5 redshift bins at 
1.0 < z < 4.0 where the VVDS AGN sample consists of 121 
objects. Since, in this redshift range, our data cover only the 
faint part of the luminosity function, these fits have been per- 
formed by combining our da ta with the LF data poin ts from the 
SDSS data release 3 (DR3) dRichards et all l2006bl) in the red- 



shift range < z < 4. The advantage of using the SDSS sample, 
rather than, for example, the 2dF sample, is that the former sam- 
ple, because of the way it is selected, probes the luminosity func- 
tion to much higher redshifts. The SDSS sample contains more 
than 15,000 spectroscopically confirmed AGN selected from an 
effective area of 1622 sq.deg. Its limiting magnitude (i < 19.1 for 
z < 3.0 and i < 20.2 for z > 3.0) is much brighter than the VVDS 
and because of this it does not sample well the A GN in the faint 
part of the luminosity function. For this reason. [Richards et al.l 
(I2006bl) fitted the SDSS data using only a single power law, 
which is meant to describe the luminosity function above the 
break luminosity. Adding the VVDS data, which instead mainly 
sample the faint end of the luminosity function, and analyzing 
the two samples together, allows us to cover the entire luminosity 
range in the common redshift range (1.0 < z < 4.0), also extend- 
ing the analysis at z < 1.0 where only SDSS data are available. 

The goodness of fit between the computed LF data points 
and the various models is then determined by the x 2 test. 

For all the analyzed models we have parameterized the lu- 
minosity function as a double power law that, expressed in lumi- 
nosity, is given by: 



<&(L,z) 



(L/L*)- a + (L/L*y 



(5) 



where <f>* L is the number of AGN per Mpc 3 , L* is the characteris- 
tic luminosity around which the slope of the luminosity function 
is changing and a and fi are the two power law indices. Equation 
|5]can be expressed in absolute magnitude^ as: 



0(M,z) = 



lQ0A(a+l)(M-M") + 1 Q0.4(/3+l )(M-M*) 



(6) 



7.1. The PLE and PDE models 



The first model that we tested is a Pure Luminosity Evolution 
(PLE) with the dependence of the characteristic luminosity de- 
scribed by a 2nd-order polynomial in redshift: 



M*(z) = AT (0) - 2.5(k lZ + k 2 z 2 ). 



(7) 



Following the finding by iRichards et ail (l2006bl) for the SDSS 
sample, we have allowed a change (flattening with redshift) of 
the bright end slope according to a linear evolution in redshift: 
a(z) = a(0) + A z. The resulting best fit parameters are listed 
in the first line of Table [2] and the resulting model fit is shown 
as a green short dashed line in Figure [9] The bright end slope a 
derived by o ur fit (q , vvds--3 . 19 at z= 2.45) is consistent with the 
one found bv lRichards etaTl (l2006bl) (a*™* = -3.1).0 

This model, as shown in Figure|9] while reproduces well the 
bright part of the LF in the entire redshift range, does not fit the 
faint part of the LF at low redshift (1.0 < z < 1.5). This appears 
to be due to the fact that, given the overall best fit normalization, 
the derived faint end slope (J3 =-1.38) is too shallow to reproduce 
the VVDS data in t his red shift range. 

IRichards etall d2005l) , working on a combined 2dF-SDSS 
(2SLAQ) sample of AGN up to z = 2.1. found that, fixing 
all of the parameters e xcept /3 and the normalization, to those 
of ICroom et alj d2004l) . the resulting faint end slope is f3 — 
-1.45 ± 0.03. This value would describe better our faint LF at 
low redshift. This trend suggests a kind of combined luminosity 
and density evolution not taken into account by the used model. 



2 <b* M = <b' L L* ■ \ln\Q-° A \ 

3 in their parameterization Ai 



-0A(a+ 1) =0.84 
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For this reason, we attempted to fit the data also including a term 
of density evolution in the form of: 



$m(z) = 3>m(0) ' 10' 



k w z+k 2 . 



(8) 



In this model the evolution of the LF is described by both 
a term of luminosity evolution, which affects M*, and a term of 
density evolution, which allows for a change in the global nor- 
malization <t>*. The derived best fit parameters of this model are 
listed in the second line of Table [2] and the model fit is shown 
as a blue long dashed line in Figure [9] This model gives a better 
X 1 with respect to the previous model, describing the entire sam- 
ple better than a simple PLE (the reduced x 2 decreases from ~ 
1.9 to ~ 1.35). However, it still does not satisfactorily reproduce 
the excess of faint objects in the redshift bin 1.0 < z < 1.5 and, 
moreover, it underestimates the faint end of the LF in the last 
redshift bin (3.0<z<4.0). 

7.2. The LDDE model 

Recently, a growing number of observations at different red- 
shifts, in soft and hard X-ray bands, have found evidences 
of a flattening of the faint end slope of the LF towards high 
redshift. This trend has been described through a luminosity- 
dependent density evolution parameterization. Such a param- 
eterization allows the redshift of the AGN density peak to 
change as a function of luminosity. This could help in explain- 
ing the excess of faint AGN found in the VVDS sample at 
1.0 < z < 1.5. Therefore, we considered a luminosity depen- 
dent densit y evolution model (LDDE), as compu t ed in the major 
X-sur veys dMivaii et al.ll2000l: lUeda et all 120031: iHasinger et all 



120051) . In particular, following Hasinger et alj ( 2005b . we as 
sumed an LDDE evolution of the form: 



<X>(M B ,z) = <D(M, 0) * e d (z, M B ) 
where: 

(1+zF 1 



<?d(z, M B ) 



along with 



Zc(M B ) 



(z < Zc) 

e d (zc)[(l+z)/(l+z c )F 2 (Z>Z C ) 



ZcolO- - 4 ^ 8 -^ (M B > M c ) 
Zc,o (M B < M c ) 



(9) 



(10) 
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where z c corresponds to the redshift at which the evolution 
changes. Note that z c is not constant but it depends on the lu- 
minosity. This dependence allows different evolutions at differ- 
ent luminosities and can indeed reproduce the differential AGN 
evolution as a function of luminosity, thus modifying the shape 
of the luminosity function as a function of redshift. We also con- 
sidered two different assumptions for pi and p2: (i) both param- 
eters constant and (ii) both linearly depending on luminosity as 
follows: 



pi (M B ) 
P 2(M B ) 



P 2 M tc 



0.4 ei (M B 
0.4e 2 (M B 



M ref ) 
M ref ) 



(12) 
(13) 



The corresponding^- 2 values for the two above cases are re- 
spectively x 2 =64.6 and;f 2 =56.8. Given the relatively small im- 
provement of the fit, we considered the addition of the two fur- 
ther parameters (ei and ej) unnecessary. The model with con- 
stant pi and p2 values is shown with a solid black line in Figure 
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Figure 10. Evolution of comoving AGN space density with red- 
shift, for different luminosity range: -22.0< M B <-20.0; -24.0< 
M B <-22.0; -26.0< M B <-24.0 and <-26.0. Dashed lines 
correspond to the redshift range in which the model has been 
extrapolated. 



[9] and the best fit parameters derived for this model are reported 
in the last line of Table [2] 

This model reproduces well the overall shape of the luminos- 
ity function over the entire redshift range, including the excess 
of faint AGN at 1 .0 < z < 1 .5. The x 1 value for the LDDE model 
is in fact the best among all the analyzed models. We found in 
fact a^- 2 of 64.6 for 67 degree of freedom and, as the reduced^ 2 
is below 1, it is acceptable^. 

The best fit value of the faint end slope, which in this model 
corresponds to the slope at z = 0, is f5 =-2. 0. This value is consis- 
tent with that derived bv lHao et al.l d2005l) studying the emission 
line luminosity function of a sample of Seyfert galaxies at very 
low redshift (0 < z < 0.15), extracted from the SDSS. They in 
fact derived a slope {S ranging from -2.07 to -2.03, depending on 
the line (Ha, [On] or [Oiii]) used to compute the nuclear lumi- 
nosity. Moreover, also the normalizations are in good agreement, 
confirming our model also in a redshift range where data are not 
available and indeed leading us to have a good confidence on the 
extrapolation of the derived model. 



8. The AGN activity as a function of redshift 

By integrating the luminosity function corresponding to our best 
fit model (i.e the LDDE model; see Tabled, we derived the co- 
moving AGN space density as a function of redshift for different 
luminosity ranges (Figure [TOb. 

The existence of a peak at z~ 2 in the space density of bright 
AGN is known since a long time, even if rarely it has been possi- 
ble to precisely locate the position of this maximum within a sin- 
gle optical survey. Figure[l0]shows that for our best fit model the 
peak of the AGN space density shifts significantly towards lower 



4 We note that the reduced x 1 °f our best fit model, which in- 
cludes also VVDS da ta, is significantly better than that obtained by 
iRichards et all d2006bl) in fitting only the SDSS DR3 data. 
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Figure 9. Filled circles correspond to our rest-frame B-band luminosity function data points, derived in the redshift bins 1.0 < z < 
1.5, 1.5 < z < 2.0, 2.0 < z < 2.5 , 2.5 < z < 3.0 and 3.0 < z < 4.0. Open circles are the data points from the SDSS Data Release 3 
(DR3) bv lRichards et~ai1(l2006bl). These data ar e shown also in two redshift bins below z = 1 . The red dot-dashed line corresponds to 
the model fit derived bv lRichards et al.l (l2006bl) only for the SDSS data. The other lines correspond to model fits derived considering 
the combination of the VVDS and SDSS samples for different evolutionary models, as listed in Tableland described in Section|7] 



redshift going to lower luminosity. The position of the maximum 
moves from z~ 2.0 for M B <-26.0 to z~ 0.65 for -22< M B <-20. 

A similar trend has recently been fo und by the analysis 
of several dee p X-ray selected samples (ICowie et all 120031 : 
iHasinger et all 120051: Ea Franca et all 12005^ To compare with 
X-ray results, by applying the same bolometric corrections 
used is Section [6] we derived the volume densities derived 
by our best fit LDDE m odel in the same luminosity ranges 
as lLa Franca et al.l d2005l) . We found that the volume density 



peaks at z ^ [0.35; 0.7; 1.1; 1.5] respectively for LogLx(2-iofev) 
= [42-43; 43-44; 44-44.5; 44.5-45]. In the same luminosity 
intervals, the values f or the redshift of the peak obtained by 
lLa Franca etafl d2005l) are z [0.5; 0.8; 1.1; 1.5], in good agree- 
ment with our result. This trend has been interpreted as evidence 
of AGN (i.e. black hole) "cosmic downsizing", similar to what 
has recently been ob served in the gal axy spheroid popul ation 
(ICimatti et all 120061) . The downsizing (ICowie et all 1 19961) is a 
term which is used to describe the phenomenon whereby lumi- 
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Table 1. Binned luminosity function estimate for Q.„,=0.3, Qa=0.7 and Ho=70 km ■ s _1 ■ Mpc" 1 . We list the values of Log cD and 
the corresponding l<x errors in five redshift ranges, as plotted with full circles in Figure|9]and in AM B =1.0 magnitude bins. We also 
list the number of AGN contributing to the luminosity function estimate in each bin 
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Table 2. Best fit models derived from the^- 2 analysis of the combined sample VVDS+SDSS-DR3 in the redshift range 0.0 < z < 4.0 
assuming a flat (Q m + Qa = 1) universe with Q„, = 0.3. 



nous activity (star formation and accretion onto black holes) ap- 
pears to be occurring predominantly in progressively lower mass 
objects (galaxies or BHs) as the redshift decreases. As such, it 
explains why the number of bright sources peaks at higher red- 
shift than the number of faint sources. 

As already said, this effect had not been seen so far in the 
analysis of optically selected samples. This can be due to the 
fact that most of the optical samples, because of their limiting 
magnitudes, do not reach luminosities where the difference in the 
location of the pe ak becomes evident. The COMBO- 17 sample 
dWolf et al.Ll2003l) . for example, even if it covers enough redshift 
range (1.2 < z < 4.8) to enclose the peak of the AGN activity, 
does not probe luminosities faint enough to find a significant 
indication for a difference between the space density peaks of 
AGN of different luminosities (see, for example, Figure 11 in 
IWolf et all d2003l) . which is analogous to our Figure [10] but in 
which only AGN brighter than M ~ -24 are shown). The VVDS 
sample, being about one magnitude fainter than the COMBO- 
17 sample and not having any bias in finding faint AGN, allows 
us to detect for the first time in an optically selected sample the 
shift of the maximum space density towards lower redshift for 
low luminosity AGN. 



9. Summary and conclusion 

In the present paper we have used th e new sample of AGN, c ol- 
lected by the VVDS and presented in lGavignaud et al] d2006l) . to 
derive the optical luminosity function of faint type-1 AGN. 

The sample consists of 130 broad line AGN (BLAGN) se- 
lected on the basis of only their spectral features, with no mor- 
phological and/or color selection biases. The absence of these 
biases is particularly important for this sample because the typ- 
ical non-thermal AGN continuum can be significantly masked 
by the emission of the host galaxy at the low intrinsic luminos- 
ity of the VVDS AGN. This makes the optical selection of the 
faint AGN candidates very difficult using the standard color and 
morphological criteria. Only spectroscopic surveys without any 
pre-selection can therefore be considered complete in this lumi- 
nosity range. 

Because of the absence of morphological and color selec- 
tion, our sample shows redder colors than those expected, for 
example, on the basis of the color track derived from the SDSS 
composite spectrum and the difference is stronger for the intrin- 
sically faintest objects. Thanks to the extended multi-wavelength 
coverage in the deep VVDS fields in which we have, in addition 
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to the optical VVDS bands, also photometric data from GALEX, 
CFHTLS, UKIDSS and SWIRE, we examined the spectral en- 
ergy distribution of each object and we fitted it with a combina- 
tion of AGN and galaxy emission, allowing also for the possi- 
bility of extinction of the AGN flux. We found that both effects 
(presence of dust and contamination from the host galaxy) are 
likely to be responsible for this reddening, even if it is not pos- 
sible to exclude that faint AGN are intrinsically redder than the 
brighter ones. 

We derived the luminosity function in the B-band f or 1 < 
z < 3.6, using the usual l/V max estimator dSchmidtl 1 19681) . 
which gives the space density contributions of individual ob- 
jects. Moreov e r, usi ng the prescriptions recently derived by 
iHopkins et al.l (120071) . we computed also the bolometric lumi- 
nosity function for our sample. This allows us to compare our 
results also with other samples selected from different bands. 

Our data, more than one magnitude fainter than previous op- 
tical surveys, allow us to constrain the faint part of the luminosity 
function up to high redshift. A comparison of our data with the 
2dF sample at low redshift (1 < z < 2.1) shows that the VVDS 
data can not be well fitted with the PLE models derived by pre- 
vious samples. Qualitatively, our data suggest the presence of 
an excess of faint objects at low redshift (1.0 < z < 1.5) with 
respect to these models. 

Recently, a growing number of observations at different red- 
shifts, in soft and hard X-ray bands, have found in fact evi- 
dences of a similar trend and they have been reproduced with a 
luminosity-dependent density evolution parameterization. Such 
a parameterization allows the redshift of the AGN density peak 
to change as a function of luminosity and explains the excess 
of faint AGN that we found atl.0<z<1.5. Indeed, by com- 
bining our faint VVDS sample with th e large sample of brigh t 
AGN extracted from the SDSS DR3 dRichards et all l2006bl) . 
we found that the evolutionary model which better represents 
the combined luminosity functions, over a wide range of red- 
shift and luminosity, is an LDDE model, similar to those derived 
from the major X-surveys. The derived faint e nd slope at z=0 is 
/3 = -2.0, consistent with the value derived by lHao et al. I (12001) 
studying the emission line luminosity function of a sample of 
Seyfert galaxies at very low redshift. 

A feature intrinsic to these LDDE models is that the comov- 
ing AGN space density shows a shift of the peak with luminos- 
ity, in the sense that more luminous AGN peak earlier in the 
history of the Universe (i.e. at higher redshift), while the density 
of low luminosity ones reaches its maximum later (i.e. at lower 
redshift). In particular, in our best fit LDDE model the peak of 
the space density ranges from z ~ 2 for Mb < -26 to z~ 0.65 
for -22 < Mb < -20. This effect had not been seen so far in the 
analysis of optically selected samples, probably because most of 
the optical samples do not sample in a complete way the faintest 
luminosities, where the difference in the location of the peak be- 
comes evident. 

Although the results here presented appear to be already ro- 
bust, the larger AGN sample we will have at the end of the still 
on-going VVDS survey (> 300 AGN), will allow a better sta- 
tistical analysis and a better estimate of the parameters of the 
evolutionary model. 
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